hep-ph/0 30202r 
IFIC/03-03 
ZU-TH02/03 



Solar Neutrino Masses and Mixing from Bilinear R-Parity 
Broken Supersymmetry: Analytical versus Numerical Results 

M. A. Dfaz\ M. Hirsch^, W. Porod^, J. C. Rom&o^'^ and J. W. F. Valle^ 

^ Facultad de Fisica, Universidad Catdlica de Chile 
Av. Vicuna Mackenna 4860, Santiago, Chile 
^ AHEP Croup, Instituto de Fisica Corpuscular - C.S.I.C./Universitat de Valencia 
Edificio de Institutos de Paterna, Apartado 22085, E-46071 Valencia, Spain 
^ Institut fiir Theoretische Physik, Universitdt Zurich, 
CH-8057 Ziirich, Switzerland 
Departamento de Fisica and CFIF, Instituto Superior Tecnico 
Av. Rovisco Pais 1, 1049-001 Lisboa, Portugal 

Abstract 

We give an analytical calculation of solar neutrino masses and mixing at one-loop 
order within bilinear R-parity breaking supersymmetry, and compare our results to 
the exact numerical calculation. Our method is based on a systematic perturba- 
tive expansion of R-parity violating vertices to leading order. We find in general 
quite good agreement between approximate and full numerical calculation, but the 
approximate expressions are much simpler to implement. Our formalism works es- 
pecially well for the case of the large mixing angle MSW solution (LMA-MSW), 
now strongly favoured by the recent KamLAND reactor neutrino data. 



1 Introduction 



Solar neutrino experiments, including the measurement of the neutral current rate for 
solar neutrinos by the SNO collaboration provide a solid evidence for solar neutrino 
conversions 2l- This has been recently confirmed by the first results from the KamLAND 

nn 

experiment using reactor (anti)-neutrinos [3|, |^ . Combining the information from reactors 
with all of the solar neutrino data leads to the best fit point ^5]: 

tan^ ^soL = 0.46, Am^o^ = 6.9 x 10"^ eV^ , (1) 

confirming that the solar neutrino mixing angle is large, but significantly non-maximal. 
The 3o" region for 6 is: 

0.29 < tan^^soL < 0.86, (2) 

based on a combination of all experimental data. However, one finds a significant reduc- 
tion of the allowed Am^Q^ range. As shown in Ref. I^i], the pre-KamLAND LMA-MSW 
region is now split into two sub- regions. At 3cr (1 dof.) one obtains 



5.1 X 10"^ eV^ < Am^Qj^ < 9.7 x 10"^ eV^ 1.2 x 10"^ eV^ < Am^^^ < 1.9 x 10"^ eV^ (3) 

Altogether, KamLAND results exclude all oscillation solutions except for the large 
mixing angle MSW solution (LMA-MSW) to the solar neutrino problem Q|. 

On the other hand, current atmospheric neutrino data require oscillations involving 
i^fi i^T 0- The most recent global analysis gives 

sin^ ^ATM = 0.5 , Aml^^ = 2.5 x 10"^ eV^ (4) 

with the 3a ranges (1 d.o.f.) 

0.3 < sin^^ATM < 0.7 (5) 

1.2 X 10"^ eV^ < Am^^^ < 4.8 x 10"^ eV^ . (6) 

These data have triggered a rush of theoretical and phenomenological papers on mod- 
els of neutrino masses and mixings, most of which introduce a large mass scale in order to 
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implement various variants of the see-saw mechanism j8|, |9|, llO| . Broken R-parity super- 
symmetry provides a theoretically interesting and phenomenologically viable alternative 
to the origin of neutrino mass and mixing Here we focus on the simplest case of su- 
persymmetry with bilinear R-parity breaking In contrast to the seesaw mechanism, 
here neutrino masses are generated at the electro-weak scale. Such low-scale schemes for 
„^ut*„ .a.se. Kave *e adva.a,e of bei,, ..able also acce.e..ato.. experiment, Q- 

neutralino 3, 3, a slepton |il2i] or a stop 

Sup ersymmetric models with explicit bilinear breaking of R-parity (BRpV) 



through the decay properties of the lightest supersymmetric particle if the LSP is a 
0, 3, 3, a slepton ^12,] or a stop QQ. 
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provide a simple and calculable framework for neutrino masses and 

n 

mixing angles in agreement with the experimental data [28||. In this model the atmospheric 
neutrino mass scale is generated at tree- level, through an effective 'low-scale" variant 
of the seesaw mechanism ll|. In contrast, the solar mass and mixings are generated 



radiatively |28|. Tree-level neutrino masses within BRpV have been treated extensively 
in the literature. 

This paper is mainly devoted to the solar neutrino masses and mixing. An accurate 
and reliable calculational method is now necessary in order to confront the model with 
the new experimental data from KamLAND and other neutrino experiments. A complete 
one-loop calculation of the neutrino-neutralino mass matrix has been given 



but is 



rather complex. On the other hand, approximations to the full 1-loop calculation which 



exist in the literature 



29| have not been tested yet against the full calculation. Especially 



in view of future experimental sensitivities we think such a "benchmark" is important. 

In this paper we give an accurate determination of neutrino mass and mixing within an 
analytical approximation and obtain formulae which can be rather simple, in some cases. 
For definiteness we will stick to the case of explicit BRpV only. This is the simplest of 
all R parity violating models. It can be considered either as a minimal three-parameter 
extension of the MSSM (with no new particles) valid up to some very high unification 
energy scale, or as the effective description of a more fundamental theory in which the 
breaking of R-parity is spontaneous Isol . Isil . Is^ . The latter implies the absence of trilinear 
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R-parity breaking parameters in the superpotential ^. 

This paper is organized as follows. In Sec. |21 we introduce the main features of the 
model and the relevant mass matrices and corresponding diagonalization matrices. In 
particular we identify the relevant Feynman graph topologies and rules, and derive ap- 
proximate formulae for the couplings relevant for the determination of radiatively induced 
solar neutrino mass scale. We give approximate formulas for the bottom quark/squark 
loop as well as for the charged scalar loop. In Sec. Elwe check the accuracy of our approxi- 
mation formulas by a comparison with a full numerical calculation, studying first the role 
of the simplest bottom-sbottom loop, then the charged scalar loop, before comparing the 
sum of the two to the full numerical result. In|3]we give simplified approximation formulas 
for the solar mass and solar mixing angle and conclude and summarize our results in El 



In this section we introduce the main features of the model and the relevant mass ma- 
trices, and develop approximate formulas, first for couplings and then for the radiative 
contributions to the neutrino masses due to the exchange of bottom and sbottom quarks, 
and due to charged scalars and charged fermion loops. 



The minimal BRpV model we are working with is characterized by the presence of three 
extra bilinear terms in the superpotential analogous to the /i term present in the MSSM 



where Wyuk includes the usual MSSM Yukawa terms, /i is the Higgsino mass term of 
the MSSM, and e, are the three new terms which violate R-Parity and lepton number. 
The smallness of may arise dynamically (the product of a Yukawa coupling times a 
singlet sneutrino vacuum expectation value) in models with spontaneous breaking of R 
parity js^. 

"'^Alternatively, such absence may arise from suitable symmetries [s^ 



2 BRpV Formalism 



2.1 BRpV Model 




(7) 
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Alternatively, the smallness of the e 



In fact any solution to the /i problem 



may arise from suitable family symmetries 



331. 



3J] potentially explains also the "ej-problem" 
In fact a common origin for the terms responsible for the explanation of the neutrino 
anomalies, and the /i term accounting for electroweak symmetry breaking can be ascribed 
to a suitable horizontal symmetry that may also predict their ratio, as in js^. 

In addition we have the corresponding soft supersymmetry breaking terms in the scalar 
potential, 

Vsoft = Ko/t + ^'^^ {-BfiH^.Ht + B.e^LtHt) (8) 

where B and the three Bi have units of mass and in Ko/* include all the usual mass 
and trilinear supersymmetry breaking terms of the MSSM. 

2.2 Rotation Matrices 

If the effective RpV parameters are smaller than the weak scale, we can work in a per- 
turbative expansion defined by <^ 1, where ^ denotes a 3 x 4 matrix given as 



~ ITli. 



2A 







C,i2 — ^TT 



2A, 





Ai 4Ao 

64 = -^K (9) 
^ 4Ao ^ ^ 

where Aq is the determinant of the 4x4 neutralino mass matrix, = g^Mi + g'^M2 
and 



Ai = fiVi + Vdei (10) 

The neutralino/neutrino mass matrix is diagonalized by a 7 x 7 rotation matrix AT 
according to 

Af*MFoX-^ = Mj!o"' (11) 
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and the eigenvectors are given by 

= Af,,^, (12) 

using the basis ip = (— iA', — iA^, if], if^, z/g, z/^, i^r). In this approximation, the rotation 
matrix can be written as 

Af*^\ (13) 

Here, is the rotation matrix that diagonahzes the 4x4 MSSM neutrahno mass matrix, 
Vy is the rotation matrix that diagonahzes the tree level neutrino 3x3 mass matrix, and 
^ij <^ 1 are the expansion parameters 0, 37 1. The terms we need are 

/ hh \ 

lfe= hh (14) 

Vai|A| aslAI az\K.\ + hh a4|A|/ 



where h = — l//i. 



gMin _ M^Vu _ M^Vd 
} (^3 — — } (^A — 



2Ao ' 2Ao ' 4Ao ' ~^ 4Ao ' 

The e parameters in eq. (fT^ are defined as = (vj^ e^, and are given by 



(15) 



ei 



62 



£3 



ee{Al + Al) - Ae(A^e^ + A^e^) 



A^e. - A^er 



(16) 



A 



Al + Al + Al 



On the other hand the chargino/charged slepton mass matrix is diagonalized with two 
different 5x5 mass matrices, 



(17) 



with the eigenvectors satisfying 



(18) 
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in the basis ip'^ = {—iX'^,H2,e~^,^~^,T^) and ip' 
Dirac fermions being 



-iX , iff , e^, /ij^, ), and with the 



\FuJ 



(19) 



To first order in the R-Parity violating parameters we have 



V 



V 



u 



R 



u 



VI* 



(20) 



where Vl* and diagonahze the charged lepton mass matrix according to V^*M^V^^ = 
^diag- For the purposes of our approximate formula, it is sufficient to take = 02x3, 



because the mixing between right-handed 



respect to by a factor of mi/MsusY 
We then have 

el = 4a 

and 



3,37 



eptons and the charginos is supressed with 



Note, that we can choose Vf* = V, 



R 



Ci2 
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02 Aj + bei 



9'^Vu 



V2A, 



2^Ah 



L3x3- 



(21) 



(22) 



where A+ is the determinant of the 2x2 chargino mass matrix. 

In the BRpV model the charged Higgs fields mix with the charged sleptons forming an 
8x8 mass matrix [28], which is diagonalized by a rotation matrix Rs±. The construction 
of R5± to first order in small (RpV) parameters is quite straightforward but lengthy. The 
interested reader can find the details in Appendix 1X1 



2.3 Approximate Couplings 

The relevant Feynman rules for the bottom-sbottom loops are, in the case of left sbottoms: 

b 

)bnb {i-l5) _j_ Qbnb (l+Ts) ' 




: \r)bnb {i--lh} i r)bnb U+75J ] 



Fl 
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with 

Oti = R],^{M.2 - \twKi) - R]2hbU:, (23) 

where ty/ = ta.n 6w- After approximating the rotation matrix A/" we find that expressions 
similar to eq. (f^!?|) with the replacement A/" — > N are valid when the neutral fermion is a 
neutralino. When the neutral fermion is a neutrino, the following expressions hold 

^ R^,h{as\A\5,, + bi,)+R%^twai\A\6,s 
Om, ~ R^,-^{ltwai-a2)\A\5,s + R%h{a3\A\5,s + bi,) (24) 

where i' = i — A label one of the neutrinos, i?^;, are the rotation matrices connecting weak 
and mass eigenstate basis for the scalar bottom quarks. In case of no intergenerational 
mixing in the squark sector R'jj^ can be parameterized by just one diagonalizing angle O^. 
The relevant Feynman rule for the charged Higgs/slepton loops is 



/ —A F/ncns (I-75) I iOcns (l+Ts) ] 




where the Of^^ji^ and O^^^ couplings are given in Appendix |Bl in Eqs. (jSElEZI)- 

After approximating the rotation matrices U, and V in the chargino sector, and A/" in 
the neutralino sector we find approximate expressions for these couplings that we will use 
below. These formulae are collected in Eqs. (jFTl EHl EDI) of Appendix IHl 

2.4 Relevant Topologies 

We now give the structure of the mass matrices relevant for the determination of solar 
neutrino masses and mixings. While in the BRpV model the atmospheric anomaly is 
explained at the tree-level, the solar neutrino masses and mixings are both generated 
radiatively. In particular, the "solar angle" has no meaning at the tree level due to the 
degeneracy of the two lightest neutrinos in this limit. 
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Diagonalizing the tree-level neutrino mass matrix first and adding then the 1-loop cor- 
rections before re-diagonalization the resulting neutrino/neutralino mass matrix has non- 
zero entries in the neutrino/neutrino, the neutrino/neutralino and in the neutralino/neutralino 
sectors. We have found that the most important part of the 1-loop neutrino masses de- 
rives from the neutrino/neutrino sector and that the 1-loop induced neutrino/neutralino 
mixing is usually subdominant. 

The relevant topologies for the one loop calculation of neutrino masses are then il- 
lustrated in FigH Here our conventions are as follows: open circles with a cross inside 




Figure 1: Topologies for neutrino self-energies in the BRpV supersymmetric model 



indicate genuine mass insertions which flip chirality. On the other hand open circles with- 
out a cross correspond to small R-Parity violating projections, indicating how much of an 
Rp-even/odd mass eigenstate is present in a given Rp-odd/even weak eigenstate. Strictly 
speaking these projections are really coupling matrices attached to the vertices, and this 
is what appears in the numerical code. However, given the smallness of Rp-violating 
effects, the "insertion-method" proves to be a rather useful tool to develop an analytical 
perturbative expansion and to acquire some simple understanding of the results. 

2.5 Bottom-sbottom loops 

The simplest contribution to the radiatively induced neutrino mass arises from loops 
involving bottom quarks and squarks is given by 3] 

%(0) = - E {o'Sloit + O'^O'^ m,5o(0, ml ml) (25) 



Bq{0, ml, 771^,) is the usual Passarino-Veltman function js^ 3^. This contribution can be 



expressed as being proportional to the difference of two Bq functions, 

AB~f'' = Bo{0, ml ml) - Bo{0, ml ml) 



as follows 



An 



Ncmb 



2srCrhlABk'' 



where we have defined 



(26) 



+ a^h (ei5j3 + ej5i^) |A| + {^l + j ^i^^j 



(27) 



(28) 



The different contributions can be understood as coming from the graphs correspond- 
ing to the first topology of Fig. ^ They have been depicted in more detail in Fig. |21 where 




as\A\5js + bij b a3\A\Si3 + bii [a2, ai)\A\Sj3 b 




Figure 2: Bottom-Sbottom diagrams for solar neutrino mass in the BRpV model 



we have adopted the following conventions: a) as before, open circles correspond to small 
R-parity violating projections, indicating how much of a weak eigenstate is present in a 
given mass eigenstate, (b) full circles correspond to R-parity conserving projections and 
(c) open circles with a cross inside indicate genuine mass insertions which flip chirality. 

The open and full circles should really appear at the vertices since the particles prop- 
agating in the loop are the mass eigenstates. We have however separated them to better 
identify the origin of the various terms. There is another set of graphs analogous to the 
previous ones which corresponds to the heavy sbottom. They are obtained from the pre- 
vious graphs making the replacement hi — > 62, si and — s^. Note that for all 
contributions to the 2x2 submatrix corresponding to the light neutrinos the divergence 



from i?o(0,m^,m? ) is canceled by the divergence from i?o(0, m^, m?^), making finite the 
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contribution from bottom-sbottom loops to this submatrix, as it should be, since the mass 
is fully "calculable" . 



2.6 Charged Scalar- Charged Fermion Loops 



Another contribution to the radiative 
scalar /charged- fermion loops, given as 

1 



n,,(o) 



induced neutrino mass comes from charged- 



(29) 



k,r 



The structure of the contribution from charged Higgs/slepton loops is substantially more 
complex than that of the bottom-sbottom loop considered above. It can be expressed as 



An,; 



167r2 



+ (i ^ j] 



(30) 



where 



Bo{0, ml m\) - Bo{0, ml ml) ; X,Y = (G^, H^, L,, L^, h, 



AB^^-^^' = 4 5o(0, ml ml) + si So(0, ml m^) - Bo{0, ml 
^^G±H±nf2 ^ cj5o(0,m2,m^±) + slBo{0,mlml^) 



4 So(0, ml ml) - 4 5o(0, ml ml 



(31) 



and 



= ^rCf {V2^7'ai|A| [gV.l.a^As - -^{ga2 + 9'a,)\A\6,. 



6, 



i3 



+ hi (be, + as\A\Si3 - cp^Vj^i^) [foe, + a^lMSjs - Vj^^^ia^A^ + fees)] } 



+ SfhlVj^[b~ej + a3|A|5,-3 - Kj3(a^A3 + ht^^)] } 



C^^'^ = spQuR, {sfhrVj^^^[gVj^^^a{A^ - j^{ga2 + g'ai)\A\6,3 
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a3\A\5js + bej ^ as\A\Sis + bei (a2, ai)|A|(5j3 ^ 



ai\A\S. 



'i3 





Figure 3: Charged scalar contributions to solar neutrino masses in BRpV model: terms 



proportional to Ai?, 



be, + a3\A\6j3 - V;^j3(«2 A3 + ^£3)] } 



a 



-SfsQHL, MK'isKjl^ A 



C^f^^ = -c^^MK;3K>iA3 
eft"-' = -cp-hr9VZsVj:,AA3 



a 



V 

Cf3— KV^^is 



1 

V2' 



gV^sO-iAs -^{90-2 + g'ai)\A\6j3 



KgeiV i^ba^ A3 



(32) 



The result of Eq. (jHn|) can be represented graphically for better understanding. The 
terms proportional to ABq^'^^ come from the graphs of Fig. IHl There is another set of 
four graphs corresponding to f2- These are found after making the replacements fi ^ 
Sf Cf, and Cf — ^ —Sf. The diagrams in the first row are the ones that are equivalent 
to those in the bottom-sbottom loop. They have as characteristic feature the presence 
of two Rp-violating insertions (open circles) in the external legs. However, in contrast to 
the quark sector, R-parity violation can also appear in the charged internal lines running 
in the loops, since it occurs in the charged fermion sector. This explains the origin of 
the second row in Fig. IHl The presence of R-parity violating insertions in the internal 
lines of the second row in Fig. El correspond to the second topology in Fig. The full 
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diagrammatic explanation of the rest of the terms appearing in Eq. ()30|) is given in detail 
in Appendix O 



3 Analytical versus Numerical results 

In this section we will check the accuracy of the approximation formulas given in Sec. 12.51 
and l2.(jl We do this by comparing the results obtained with their use with a full numerical 
calculation of the one-loop contributions to the neutrino mass, whose details can be found 



in Ref. 



2^ 



As will be explained in more detail below, the relative importance of the various loops 
depends on the - currently unknown - supersymmetric parameters. In order to reduce 
the number of free parameters in the following we will adopt the minimal constrained 
supergravity (mSUGRA) version of the MSSM. As a rule of thumb it can be said that 
the bottom-sbottom loop usually gives the main contribution to the neutrino mass matrix 
when the neutralino is the LSP. On the other hand, if the scalar tau is the LSP, both 
bottom-sbottom and charged scalar loops are of approximately comparable magnitudes. 

We have therefore constructed two different random scans over SUSY parameter space. 
Both sets start with the following rather generous parameter ranges: M2 from [0,1.2] TeV, 

from [0,2.5] TeV, mo in the range [0,1.0] TeV, Ao/mo and Bo/niQ [-3,3] and tan/? 
[2.5,10]. All randomly generated points were subsequently tested for consistency with the 
minimization (tadpole) conditions of the Higgs potential, as well as for phenomenological 
constraints from supersymmetric particle searches. We then selected points in which a) 
the lightest neutralino is the LSP (called set "Ntrl" in the following) or b) at least one of 
the charged sleptons was the LSP (called set "Stau" in the following). Note that in the 
Stau set niQ « M2 and large fi values are strongly preferred. 

R-parity violating parameters are chosen in such a way that neutrino oscillation data 
are reproduced approximately. As discussed in the introduction, atmospheric neutrino 
experiments require a near-to- maximal atmospheric mixing angle ^atmj with Am^^^j in the 
range given in Eq. (0). On the other hand reactor data constrain the electron-neutrino 
component in the third mass eigenstate to be small. And, finally, in combination with 
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solar neutrino data, the KamLAND data require a 6'goL in the range given in Eq. Q with 
ArrigQj^ as given in Eq. 0. The latter ranges belong to the LMA-MSW region indicated 
by a solar-only global analysis of neutrino data given in Ref. j^. For completeness we 
also include the (pre-KamLAND) LOW and VAC-type solutions of the solar neutrino 
anomaly. In the following we will first discuss the bottom-sbottom and the charged 
scalar loops separately, before considering a calculation taking into account both loops in 
comparison to the full calculation. 

3.1 Bottom-sbottom loop 

In Fig. 0] we show the ratio of the approximate-over-exact solar neutrino mass parameter 
^Appr Inexact ygj^g^g Amg^^^ for the case in which only the bottom-sbottom loop is taken 
into account, both in the approximate and in the exact calculation. The horizontal bands 
indicate attainable neutrino mass values when the parameters are scanned as indicated 
previously. As can be seen from the figure the approximate formula works quite well 
for points in both Ntrl and Stau sets, as long as the neutrino masses fall in the LMA- 
MSW range indicated by the right vertical bands. Note that the LMA-MSW and LOW 
bands indicated in the figure correspond to the full analysis of solar data only, presented 
in Ref. 1^. The recent KamLAND reactor neutrino data rule out the LOW solution 
and restricts the LMA-MSW to somewhat narrower ranges indicated in Eq. ©. One 
finds that the mass values inferred from our present analytical approximation are always 
within 10 % or less of the exact numerical calculation of the bottom-sbottom loop. Larger 
deviations show up only in the Ntrl set, for very small neutrino masses, which we trace 
to the neglection of the 1-loop neutrino/neutralino mixing terms in our approximate 
treatment. Although not strictly ruled out by a solar-only global neutrino data analysis P| , 
these LOW and VAC-type solutions are now strongly disfavored by the latest KamLAND 
reactor neutrino data. 
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Figure 4: Ratio m^^P'^ / ml^"'^* versus AmgQ^ iii for the sets Ntrl (left) and Stau (right), 
for a calculation involving only the bottom-sbottom loop. The vertical grid lines indicate 
the 90 % c.l. regions for the LOW and LMA solutions to the solar neutrino problem. 
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Figure 5: {m^^^^ / ml^"''^^) versus Am^Q^ [eV^^] for the sets Ntrl (left) and Stau (right), for 
a calculation involving only the charged scalar loop. 
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Figure 6: {m{^^'P'^ /rrfj^"'^^) versus ArngQ^ [e^^] for the set Ntrl (left) and the set Stau 
(right). m^^P^ is the sum of the bottom-sbottom and charged scalar loops, while m'^^^^ 
is the numerical result for all loops. In case of LMA the approximation works always 
better than 10 %. For the LOW solution the typical error is of the order of 10 %, while 
in extreme cases errors up to 25 % can be found. 

3.2 Charged scalar loop 

In Fig. El we show the ratio of the approximate-over-exact solar neutrino mass parameter 
^Appr Inexact versus ArngQ^, for a calculation which takes into account only the 

charged scalar loop in both the approximate and the exact calculation. As can be seen 
from the figure the approximate formula is accurate for all points in the LMA-MSW 
region, indicated by the right vertical bands 2], both for the Stau and for the Ntrl sets. 
The only case where our analytic results gives a poorer approximation (to better than a 
factor-of-2) of the full numerical result is for the Ntrl set, when the neutrino mass falls 
in the LOW or VAC ranges, now strongly disfavored by the KamLAND results. We have 
checked numerically that for these very small neutrino masses all terms in eq. (jHU)) are 
of approximately equal importance and there are significant cancellations among terms, 
which leads to a less reliable final result. 
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3.3 Comparison with full calculation 

In supersymmetric models with mSUGRA-like boundary conditions the bottom-sbottom 
and the charged scalar loop usually give the most important contribution to the neutrino 
mass matrix. This is demonstrated in Fig. |B1 (left) for the set Ntrl and in Fig. IHl (right) 
for the set Stau. In both figures we show the ratio of the approximate-over-exact solar 
neutrino mass parameter m^^'^'^ /rrf^"''^ versus Am^Q^ in eV^, where m^^^^ is the approx- 
imate loop calculation involving the bottom-sbottom and the charged scalar loop, while 
jy^exact jg ^Yi^ exact numerical computation taking into account all loops. 

In the region of Aml^^ appropriate for the currently preferred LMA-MSW solution 
to the solar neutrino problem one finds that the approximate calculation reproduces the 
exact result better than 10 %. Only in the set Ntrl one finds larger deviations, up to 25 % 
in extreme cases, when Am^Q^ lies in the LOW region, strongly disfavored by KamLAND. 
This is due to the larger errors in the bottom-sbottom calculation in this set for small 
neutrino masses as discussed above. 

4 Simplified approximation formulas 
4.1 The solar mass 

First we note that for nearly all points in our random sets we find that m,y^ <^ m,^^. In 
other words, bilinear R-parity breaking favors a hierarchical neutrino spectrum. Moreover, 
we have found numerically that the terms proportional to x ij in the self energies 
in Eq. (f77j) give the most important contribution to m^^ in the bottom-sbottom loop 
calculation in most points of our sets. If these terms are dominant one can find a very 
simple approximation for the bottom-sbottom loop contribution to rrii^^- It is given by 

m., ^ sin(2^^)m,Ai?^-i (33) 

We have checked numerically that Eq. (jHHj) reproduces the result of the full approxi- 
mative formula to high accuracy if mj^j < 0.3m^^ . Note also that Eq. (|33|) holds only if 
the 1-loop contributions to the neutrino mass matrix are smaller than the tree-level one. 
This condition requires that |e|^/|A| < 1 approximately, i.e. the bilinear parameters 
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must be suppressed wth respect to Note that sueh a suppression eould. in prineiple. 

be motivated by suitable flavour symmetries |33l ]. 

Due to the more complicated structure of the charged scalar loop it is not possible to 
give a simple equation for rrii^^ similar to Eq. ()33p for the bottom-sbottom loop. However, 
for mjyj larger than (few) xlO~^ we have found that the most important contributions to 
the charged scalar loop are the terms proportional to ABq^'^^ , AB^^'^^ and AB^^'^'^ in 
Eq. (jHUj) . We note in passing that Eq. (jHSl), with appropriate replacements, allows us 
to estimate the typical contributions to the charged scalar loop within a factor of ~ 3. 
However, such an estimate will be biased toward too small (large) m^^ for scalar tau 
(neutralino) LSPs. 




Figure 7: {m^^^"^ / rrfj^"''^) versus Amg^^^ [eV^^] for the set Ntrl (left) and Stau (right). Shown 
is the result of the simplified approximation formula in Eq. (j33|) for the sbottom-bottom 
loop and taking into account only coefficients Ch±t2i C*i?±fi and Cf,^f-^ in the charged scalar 
loop. 

In Fig. (|7I) we show a comparison of our simplified approximation formula, including 
the simple form of the sbottom-bottom loop and the three most important coefficients for 
the charged scalar loop, as discussed above, to the full numerical calculation including all 
loops. As one can see, even the simplified version of our formula works surprisingly well in 
the LMA-MSW regime, although the agreement with the full calculation is now less good 
for the LOW region, as could have been expected from the results discussed previously. 
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\ 



+ 



(34) 



4.2 The solar mixing angle 

In the basis where the tree-level neutrino mass matrix is diagonal the mass matrix at 
one-loop level can be written as 

cieiei cieie2 cieies 
Ci^2^i Cie2e'2 £16263 
016361 016362 Co|Ap + 016363 
where the q were defined before in Eq. (fT^. Coefficients cq and ci contain couplings 
and supersymmetric masses. Since they cancel in the final expression for the angle their 
exact definition is not necessary in the following. Dots stand for other terms which we 
will assume to be less important in the following, see the discussion at the end of this 
subsection. This matrix can be diagonalized approximately under the condition 

X = < 1 (35) 

ColAP 

i.e. if the 1-loop contribution to the neutrino mass matrix is smaller than the tree-level 
contribution, as also discussed above for Eq. (jH^ . Then 
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(36) 



We now calculate the eigenvalues and eigenvectors of this matrix as series expansions in 
the small x parameter. For the eigenvalues we get 





mi 
m2 



X Cq- 



0{x^) = ci {el + et)+0{x^ 



(37) 



7713 = Co|Ap + Ci?2 + C(x2) 

and for the first two eigenvalues (the third can also be easily obtained but it will not be 
necessary for the discussion of the solar mixing angle), 

62 



ei 



62 



61 



6? 



61^61 + 61' \ei + e 

(,62,1,62,2,62,3) 
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, 



(3J 



where up to O(x^) we have, 



62,1 



62,2 



62,3 



6163 



el(ef + el)^2 



1 eie3Ve3(?i + ^) 



6263 



1 6263 Je|(ef + el) 



(39) 



a; H + 



l^r 1^1 
Knowing the eigenvectors we can write down the rotation matrix that diagonahzes rhy, 



V^niyVy = diag(mi, 777-2, m3) 



(40) 



where 



61,1 61,2 61,3 



62,1 62,2 62,3 
y 63,1 63,2 63,3 J 

The neutrino mixing matrix is then given by 



(41) 



u = ivjvj 



(42) 



Using the fact that Ugs has to be small one can get the following expression for the solar 
mixing angle: 



tan^ 6'soL 



TP 

TP 



Now using the Eqs. (jUJ, (jH^ and substituting in Eq. 
expression for the solar mixing angle. 



2 

tan ^soL = :3T 
62 



(43) 

we obtain the very simple 

(44) 



This formula is a very good approximation if the one-loop matrix has the structure x ej, 
as is the case of the bottom-sbottom loop (and, to a lesser extend also for the charged 
scalar loop, which has one coefficient with the same index structure), and if 777^3 ^ m^^ . 
This is illustrated in Fig. |H| 
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In the left panel we show a calculation comparing for all points in the set Ntrl the 
approximate to the exact solar angle, while the right panel shows a subset of points using 
the cut sin(26'^)Ai?g^'^^ > 0.02. Note that this cut is designed such as to prefer points in 
which there is a sizeable contribution to the full 1-loop neutrino mass due to the bottom- 
sbottom loop. For points in which the charged scalar loop dominates eq. (j44|) gives only 
a factor-of-two estimate of the true solar angle. 

Note finally that eq. will fail completely, if = A,- and = e ^, s ince then 



e| = 0, see Eq. |T6|l . This is the origin of the "sign condition" discussed in 



< o 
<x> 

CM 

c 




2^Exact 

tan e^Q^ 



<x> 



a 



< o 
<x> 
a 

03 




0.2 0.4 0.6 0.8 1 



2„Exact 



Figure 8: (tan2^soL''''''7tan2^ 

^^^exac ^ ygj^g^g tan^ ^soL*^^"^^- On the left panel the darker 
region contains over 90% of the points in our sample. In the right panel the points in the 
region shown satisfy the cut sin(2^^jj)Ai?g^'^^ > 0.02 . 
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5 Discussion and conclusions 

We have presented an approximate calculation of the neutrino mass matrix at one-loop 
in supersymmetry with bilinearly broken R-parity. The method is based on a systematic 
perturbative expansion of R-parity violating vertices to leading order. We have identified 
the bottom-sbottom and the charged scalar loop as the most important ones, at least in 
supersymmetric models with mSUGRA-like boundary conditions. Taking into account 
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only these loops, we have given explicit formulas and discussed their validity as well 
as the accuracy with which they describe solar neutrino mass and mixing parameters. 
This was done by comparing our analytical results to the exact numerical calculation. 
We have found that for the case of the large mixing MSW solution our formulas - even 
within the simplified form Eq. (j33|) and Eq. (j44p - yield good agreement with the full 
numerical calculation, but are much simpler to implement than the full numerical one-loop 
calculation. The only solar neutrino "solutions" for which our analytical approximation 
is less accurate are those that are now ruled out by the recent reactor neutrino data from 
KamLAND. 

Let us finally discuss some possible caveats to the success of our approximate treat- 
ment. One is the assumption that supersymmetry breaking mass terms are flavour diag- 
onal, which we have adopted, motivated by constraints from flavour changing processes. 
Although such terms could be included into our approximate treatment, we have not 
done so, mainly due to the fact that the resulting formulas would be much more compli- 
cated and, therefore, of very limited practical use. A second concern is that our sample 
points were all generated using mSugra assumptions for the soft breaking masses. Clearly 
there are other possibilities to break supersymmetry and even though we expect that 
the bottom-sbottom loop and the charged scalar loop will still be well described by our 
approximation formulas, other loops, which we didn't take into account, might be more 
important than what we have found in our data sets. 
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A Rotation Matrices 



G 



(45) 



In the basis {H^,H^,e']^,ii'^,T^,e~^,ii~^,T^), one can write, to first order in R-parity 
violating parameters, the Goldstone rotation matrix as 

C/3 —Sff Vi/v V2/V vs/v 0" 

sp 0/3 

-c^vi/v spvi/v 1 

-C13V2/V S13V2/V 1 

-c^vs/v spvs/v 1 

1 

1 

000 1. 

where v'^ = v'^ + to this order, and tan/3 = Vu/vd, as usual. We have also used the 
shorthand notation 0/3(5/3) = cos/3(sin/9). 

Neglecting the electron and muon Yukawa couplings, the rotation that diagonalizes 
the sleptons at tree level is given by (in the same basis as above) 
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(46) 



After the rotations RfR-c are performed, the charged scalar mass matrix is diago- 
nalized up to small R-parity violating entries. In the approximation where there is no 
intergenerational mixing and /i^ ~ /ig ~ 0, these are 
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where 



^HL3 



HRi — ^HRi, {i — 1, 2) 
HRs = —Sf^HLs + CfXnRs 



with 



and 



^HLi — SpXuLi + CpXdLi XhR3 — SfjXuR^ + CpXdR^, (i— 1)3) 



XuLi = VdVi - ll€i - ^hlvdViSis , XdLi 
XuR3 = -^hr{ArV3 + CsVu) , XdR^ = -^KinV^ + €3^^) . 

These mixings are removed with the rotation matrix Rx given by 
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in the small mixing approximation sin © ~ 0. Note that here we have defined 



m]j± - m| 
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Putting everything together we get the final form of the charged scalar diagonalization 
matrix RxRfR-G which can be expressed as 



C/3 




vi/v 


V2/V 
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-cpOnLi + S/3 V 
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Cf . 



where we have defined, 

^HLs — Cf^HLa — SfOnRs , ^HRs — Sf^HLs + Cf^HRs ■ (54) 
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B Charged Higgs/slepton couplings 

The couplings of the five (generahzed to include also the three charged leptons) charginos 
to the eight charged scalars (including Higgs bosons and sleptons of both chiralities) and 
seven neutralinos (generalized to include also the three neutrinos) are given by 2^ 

OZ]k = RlthrM,,V,,-RS{^-^M,^V,, + ^M,,^^^ (55) 

where i labels the charginos, j labels neutralinos, and k labels the charged scalars, re- 
spectively. For the the right-handed couplings the corresponding couplings are given by 



g g 

A/'j2Wi4 + -j=AfjiUi4 - gAfjeUil 



(56) 



+RisK {MjiUi2 - ATjsW.s) 

After approximating the rotation matrices W, and V in the chargino sector, and M in 
the neutralino sector we find the expressions given in Eqs. (j57|) - (pTn|) . Note that we have 
divided them into cases where the charged fermion is a lepton or a chargino. For the left 
couplings when the charged fermion is a chargino we have, 

^ |A|5,3 (57) 

where V is the reduced 2x2 chargino diagonalization matrix of the MSSM, and = 1, 2. 
If the charged fermion is a lepton we have 



Ocns 
Lijk 



Rk2 



Ocns 
Lijk 



+ 



Rk6^^l + Rkt^i2 + Rks^i^ V2g'ai\A\Sjs 



5±, 



(5J 



For the right-handed couplings when the charged fermion is a chargino we get 



Ocns 
Rijk 



R[ 



5± 
kl 



-^{90.2 + g'ai)\K\5j:iUi>2 + giptj + a3|A|(5j3)[/i/i 
RltgyJ.iUi'i - Rlt9Vj^j2U,i - RSgVZsU.i + RltKVj^,,U,2 (59) 
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where U is the second 2x2 chargino rotation matrix of the MSSM. Finally, if the charged 
fermion is a lepton one has 



^Rijk 



-R 



-R 



-R 



s± 

fc3 

S± 
■kA 

■fe5 



75(^/02 + g' ai)\K\5j^5ii - gVjj^a^K 
73(5'«2 + g'ai)\A\5j^5i2 - gYj^j^^i^i 
^{ga2 + g'ai)\A\5j-i5i3 - gVj^j^a^Ai 



(60) 



'i3 
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C Charged Scalar— Charged Fermion Loops 



There are nine different terms contributing to the Charged Scalar-Charged Fermion loop, 
as it was shown in Eq. ()30p . All these terms give a finite contribution to the 2x2 submatrix 
corresponding to the light neutrinos. In this Appendix we will explain with graphs the 
origin of the different terms. The conventions used were explained in section] 



C.l AB^'^' 

The terms proportional to AB^^^^ come from the graphs of Fig. El as explained in sec- 
tion |2J 



C.2 ABI^ and AB^^ 

Now consider the terms proportional to ABq^'^^ and ABq^^^ in Eq. Of these terms, 
the ones which are related to the charged Higgs mixing with staus, can be understood 
as coming from the four graphs of Fig. El Associated to these charged Higgs graphs are 
those related to the fi mixing with charged Higgs. These are given in Fig. E3 




O—'^ \hr .or I 



a3\A\Sj3 + be 






a^Aa + kg- 



Figure 9: if^ contribution to ABq^'^^ 



There is another set of four graphs corresponding to T2 that are obtained from those in 
Fig. ^Iby replacing fi f2, Sf — * Cf and Cf —Sf. These three groups of four graphs. 
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[a2,ai 



IMS. 



'j3 



r 




a3\A\Sj3 + bej ^ 



Figure 10: n contributions to AE^*^^ and AB^^^^^^ 



when combined, form a set which is ultraviolet finite and account for the terms in Eq. (|3U|) 
proportional to A5^^^^ and A5^'^^^ 



C.3 AB^^^' and AB^^^' 




Figure 11: a) contribution to ABq^^^; b) cl contribution to ABq^^^ and ABq^^^ 



We now turn our attention to the terms proportional to ABq and Bq ^'^ which 
are related to the mixing between charged Higgs with selectrons and smuons. The terms 
proportional to ABq come from the diagrams of Fig. ^2 The terms proportional to 
ABq^^^ are easily obtained from these by replacing the corresponding slepton lines and 
couplings. Notice that in Fig. [TT]b) there is a contribution proportional tovi/v that does 
not belong to this term. We will show below that it will contribute to the ABq ^ term. 
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C.4 /\bT^' and /\bT^' 




Figure 12: contribution to AB^^^^ 

The graphs contributing to the ABq and ABq terms are related to those of 
Fig. ^2 They are given by Fig. IT^ and by the term proportional to vi/v in Fig. ITD b). 
for the case of the selectron. The terms proportional to ABq are easily obtained from 
these by replacing the corresponding slepton lines and couplings. 



C.5 AB^ 

We now consider a more complicated term, the one proportional to ABq This term 

gives a finite ultraviolet contribution and comes from the diagrams of Fig. together 
with the parts of the diagrams of Fig. ^Jthat are proportional to v^/v. Corresponding 
to the diagrams in Fig. proportional to v^/v, there is another set with fi and f2 
interchanged in the usual way. 



V 




Vi 




V, i3 



(a2,ai)|A|Jj3 



Figure 13: contribution to AB^^^^^^ 



C.6 AB^^^^^'^' 

Let us consider finally the last term in Eq. (j3U|) . the one proportional to ABq'^^^'^^'^^ . 
This term gives an ultraviolet finite contribution and comes from four diagrams. The first 
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two are those represented in Fig. corresponding to a iJ^ and ti propagating in the 
loop. The other two are obtained from these with the replacements, 



, S(3 Cf3 



(61) 



6ei 




d. 



Figure 14: and Ti contributions to ABq^^^'^^'^^ 
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